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Abstract. We discuss evidence in local, present-day clusters ol galaxies (from the ENACS sur- 
vey) about the way in which those clusters have evolved and about the evolutionary relationships 
between the galaxies of different morphological types in them. This evidence is complementary 
to that obtained from the study of clusters at intermediate and high redshifts. We argue that 
the spatial distribution and the kinematics of the various types of galaxies in and outside sub- 
structures support the following picture. 

The elliptical and SO galaxies have been around for a long time and have obtained an isotropic 
velocity distribution. The spatial distribution and kinematics of the early spirals are consistent 
with the idea that many of their kind have transformed into an SO, but that they have survived, 
most likely because of their velocities. The distribution and kinematics of the late spirals are 
consistent with a picture in which they have been accreted fairly recently. They have mildly radial 
orbits and hardly populate the central regions, most likely because they suffer tidal disruption. 
Finally, the distribution and kinematics of the galaxies in substructures, when taken at face 
value, imply tangential velocity anisotropy for these galaxies, but this result may be (partly) 
due to the procedure by which these galaxies are selected. A first attempt to take the effects of 
selection into account shows that isotropic (or even mildly radial) orbits of subcluster galaxies 
cannot be excluded. 



1. Introduction and summary of the data sample 

By observing the precursors of present-day clusters, it has become possible to study in a 
direct manner the evolution of the clusters themselves, and the evolutionary relationships 
between the various types of galaxies in them. Images of rich clusters at redshifts z « 1 
show a high merger fraction, which is taken to be direct evidence for the formation 
of elliptical galaxies by merging. Yet, it is unlikely that all early- type galaxies formed 
around or before that epoch (van Dokkum & Franx 2001). In particular, the variation 
with redshift of the composition of the galaxy population in clusters suggests that SO 
galaxies have come into being fairly recently, say, since z ~ 0.5, apparently at the expense 
of the spirals (Dressier et al. 1997; Fasano et al. 2000). The processes by which spirals 
transform into SO galaxies have been studied numerically (e.g. Moore et al. 1998; Abadi 
et al. 1999). It appears that impulsive encounters can indeed transform a spiral into an 
SO galaxy, by stripping a small fraction of the stellar disk, and by heating up the disk. 

We have used the data from the ENACS to investigate if the effects of these processes 
are still detectable in present-day clusters. Our analysis is based on an 'ensemble' clus- 
ter of ~ 3000 galaxies built by combining 59 nearby (0.035 ^ z ^ 0.098) clusters from 
the ENACS data-set (Katgert et al. 1996, 1998), after suitable scaling of the galaxies 
projected clustercentric distances and velocities (Biviano et al. 2002). For most of these 
galaxies also spectral- type or morphological data are available (Thomas 2002). We con- 
sider five populations of cluster galaxies, selected on the basis of their different projected 
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phase-space distributions (Biviano et al. 2002). These populations are: (i) galaxies in 
subclusters (herafter referred to as the 'Subs' class), and, among the galaxies outside 
subclusters: (ii) the brightest ellipticals (with Mr ^ —22 + 51og/i, the 'Bright' class), 
(hi) the other ellipticals together with the E/SO and the SO galaxies (the 'Early' class), 
(iv) the early spirals (Sa-Sb; the 'S-early' class), and (v) the late spirals and irregulars 
(Sbc-Ir) together with the emission-line galaxies (except those with early morphology), 
hereafter globally referred to as the 'S-late' class galaxies. 



The substructure distinction is based on a refined version of the Dressier & Shectman (1988) 



criterion. Galaxies in cold and/or moving substructures were identified by the velocity dis- 
persion and/or average velocity of their surroundings, which are smaller than the global 
value, or different from the cluster mean, respectively. With the significance threshold 
that we used, about 25 % of all galaxies are in substructures. The fraction of galaxies 
in substructures shows a significant apparent decrease towards the cluster center. This 
is probably a real effect, but there is likely to be a systematic contribution from the 
selection algorithm (Biviano et al. 2002). 



2. Clues from the spatial distributions 

The projected distributions of the 5 classes are shown in the left-hand panel of Fig. 1, 
in the form of the smooth estimate of I(R) (using the LOWESS technique, see Gebhardt 
& Fisher 1995; Katgert et al. 2004; Biviano & Katgert 2004). Note that our cluster 
sample has a median value r2oo ~ 1.2 ft, -1 Mpc. The brightest ellipticals (Bright-class 
galaxies) are most centrally concentrated. The Early-class shows a flatter slope, however 
with a clear steepening beyond « 0.3 7-200 ■ The S-early- and S-late- classes both avoid 
the central region within w 0.1 T200i but in addition the S-late-class shows an appreciable 
flattening of I{R) below 0.3r2oo- The Abel deprojection of I(R) which yields the 3-D 
number density profiles shows that both the S-early- and S-late- classes are absent from 
the central region. Apparently, galaxies of the S-early- and S-late- classes get destroyed or 
transformed within « 0.3r2oo, as it is unlikely that their phase-space distributions would 
prohibit them from entering the central regions. Galaxies in substructures do not avoid 
the central regions, but are less centrally concentrated than galaxies of the Early-class. 

Information on possible transformation relations between the galaxies in the differ- 
ent classes may be gleaned from the morphology-radius (MR) and, in particular, the 
morphology-density (MD) relation (e.g. Dressier 1980). For the EN ACS clusters those 
have been discussed by Thomas & Katgert (2004) who considered ellipticals (including 
the brightest) and SO galaxies separately, in addition to the S-early- and S-late-classes. 
The MR-relations show the dependence on radius of the composition in terms of mor- 
phology, and are fully consistent with the I(R)'s in the upper panel of Fig. 1. The 
MD-relation gives the distribution of local projected density around galaxies of the dif- 
ferent classes. Thomas fc Katgert (2004) first used the original definition of projected 
density, i.e. that based on the 10 nearest neigbours, or £10. However, £10 contains a very 
noticeable print-through of the average relation between projected density and radius. As 
a result, the £10 MD-relation is highly correlated with the MR-relation (and thus, with 
the I(R)'s). The ellipticals prefer the highest densities, and the remaining three classes 
have indistinguishable £10 distributions shifted towards lower densities (left-hand panel 
of Fig. 2). From this perspective, the SO galaxies and ellipticals behave differently; part 
of this difference is due to the brightest ellipticals being located at the clusters centres. 



Thomas & Katgert (2004) also used £1, i.e. the projected density based on the distance 



to the nearest neighbour. As expected, £1 is much noisier than £10, but at the same time 
it is much less correlated with projected radius, and therefore a less biased measure of 
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Figure 1. The projected distributions (left) and the velocity dispersion profiles (right) of the 
five galaxy classes. Dashed lines indicate 1-cr confidence levels. Dots (with 1-a error bars) in the 
right-hand panel indicate binned estimates of the velocity dispersion profiles. Units along the 
y-axis in the left-hand panel are arbitrary. The y-axis in the right-hand panel is in units of the 
cluster global line-of-sight velocity dispersion (on average, « 700 km s _1 for the 59 clusters that 
compose the ensemble cluster). 



local projected density than Sio. That is probably the reason why the Si-distributions 
of SO and S-early galaxies are significantly different, with SO galaxies preferring higher 
densities than the S-early galaxies (right-hand panel of Fig. 2). 
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3. Clues from the kinematics 

3.1. Isotropic orbit distributions? 

For a study of the kinematics of the various galaxy classes, we have used the observed 
projected distributions I(R) and velocity dispersion profiles (vdp's, hereafter). The vdp's 
of the five different classes are quite different (see right-hand panel of Fig. 1). The Bright- 
class vdp is quite 'cold', not only in the cluster centre. The Early-class vdp is nearly flat. 
The vdp's of the S-early and the S-late classes are rather similar, showing a steep decline 
from the centre out to ~ 0.3r2oo, followed by a nearly flat vdp. However, the S-late vdp 
is 'hotter' than the S-early vdp at all radii. Finally, the vdp of the Subs-class is cold and 
flat, even flatter than the Early-class vdp. 

In order to determine the kinematics of the different galaxy classes we first need to 
determine the gravitational potential in which the galaxies move. We have therefore first 
derived the mass-profile M(< r), using the data for the Early-class galaxies, assuming 
that those have an isotropic distribution of orbits. That assumption is justified by the 



con- 



shape of their line-of-sight velocity distribution, from which |Katgert et al. (2004) 
eluded that, for the Early-class 0.8 < j3' < 1.05, with 

/?' = ^ r Rfl\ (3.1) 

where v^, v\ are the mean squared components of the radial and tangential velocity, 
respectively (i.e. isotropic orbits correspond to (3'{r) = 1). Therefore the cluster mass 
profile, M(< r), follows from the isotropic Jeans equation (Binney & Tremaine 1987) 
applied to the Early-class. The resulting M(< r) is rather similar to a NFW (Navarro et 
al. 1997) profile with a concentration parameter c = 4.0^'g (Katgert et al. 2004). 

Next, this cluster M(< r) was used to search for solutions of dynamical equilibrium 
for the other four classes (see §0. Given the cluster M(< r) and the observed number 
density profile of a given class of cluster galaxies, one can solve the inverse Jeans equation 
(van der Marel 1994) for a given (3'(r). The resulting vdp is then compared, after the 
usual Abel projection (Binney & Tremaine 1987), with the observed line-of-sight vdp, in a 
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X 2 sense. If the two vdp's are statistically indistinguishable, the f3'(r) profile is considered 
acceptable. Instead of considering complicated f3'(r) profiles, we tested whether the data 
were consistent with an isotropic orbit distribution for the other galaxy classes, i.e. if 
P'(r) = 1 is acceptable. This appears to be the case for galaxies of the S-early class, but 
not for those of the Bright, S-late, and Subs classes. 

3.2. The orbital anisotropy profiles 

Since the assumption of (nearly) isotropic orbits is not valid for all except the Early and, 
possibly, the S-early classes we have investigated what the data tell us about the orbital 
anisotropy. Instead of trying constant anisotropy (van der Marel et al. 2000; Lokas & 
Mamon 2003), or anisotropy models (Carlberg et al. 1997), we determine (3'(r) directly 



from t he data, using the procedu re that Solanes & Salvador-Sole (1990) developed from 
that of |Binney fc Mamon (1982)] The f3'{r) solutions are checked by inverting the Jeans 
equation again (van der Marel 1994), thereby deriving a vdp that can be compared with 
the observed one. The derived vdp's always agree very well with the observed vdp's. 
We also applied this check to the sample of Early-class galaxies, and indeed recovered 
13' (r) = 1 to within ±15%. 

Uncertainties on the /3'-profiles were estimated by dividing each galaxy sample into 
four subsamples, each containing a different half of all the galaxies in the original sample. 
We then applied the procedure of |Solanes fc Salvador-Sole (1990)| to each subsample. 
The rms of the four resulting /3'-profiles is an approximate, but probably conservative, 
estimate of the uncertainty in the full sample (3'{r). 

Physically acceptable solutions are found for the Early and three of the four other 
populations of cluster galaxies. No acceptable solution is found for the Bright class. This 
is most likely the consequence of a failure of the collisionless Jeans equation in describing 
the dynamics of these galaxies, which is probably affected by dissipative processes, such 
as mergers and dynamical friction (see Fusco-Femiano & Menci 1998 for a collisional 
solution) . 

The /3'-profile obtained for the S-early class increases rapidly from /3' ~ 1 near the 
centre, to (3' ~ 1.8 at r/r2oo — 0.45, and then smoothly decreases to ~ 1 again at 
larger radii (see Fig. 3, top panel). Taken at face value, this solution implies mild radial 
anisotropy just outside the cluster centre. The partial radial anisotropy at r ~ 0.45 r2oo 
could result from a natural selection effect. I.e., among S-early galaxies, those reaching 
the high-density inner cluster region can avoid transforming into SO's only if they move 
sufficiently fast. We should not overinterpret the data, though, since the isotropic solution 
is acceptable for galaxies of the S-early class, and the /3'-profile is in fact consistent with 
(3'(r) = 1 . 

The (3'{r) profile of the S-late galaxies is not different from unity out to ~ 0.7r2oo, 
but then increases linearly with radius reaching a radial anisotropy ~ 1.8 at r ~ 
1-5 ^200 (see Fig. 3, mid panel). This /^'-profile is remarkably similar to those obtained 
for dark matter particles in numerical simulations (Ghigna et al. 1998; Diaferio 1999). 
It is therefore tempting to conclude, by analogy, that S-late galaxies, like dark matter 
particles in numerical simulations, still retain memory of the process of (mostly radial) 
gravitational infall along the filaments connecting to the cluster, and are therefore recent 
arrivals into the cluster. 

Finally, the data for the galaxies of the Subs class imply orbits with substantial tan- 
gential anisotropy at all radii (see Fig. 3, solid line in the bottom panel). However, two 
systematic effect must be accounted for. The first one is the internal velocity dispersion 
of subclusters. If we subtract an assumed internal velocity dispersion of 250 km s _1 
from the observed vdp of Subs galaxies, the orbital tangential anisotropy becomes even 
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Figure 3. Top panel: The velocity anisotropy profile, f3'(r) = (« 2 (r)/t)f (r)) 1 ^ 2 , of the S-early 
class. The hatched region represents the approximate 1-cr confidence interval. Mid panel: same 
as top panel, but for the S-late class. Bottom panel: same as top panel, but for the Subs 
class. The dashed line is the /^'-profile obtained after subtracting in quadrature the internal 
velocity dispersion of the subclusters from the observed vdp. The dash-dotted line is the /^'-profile 
obtained after correcting the observed number density profile and vdp of the Subs class for the 
systematic selection effects, according to a plausible, but by no means unique, model. 

stronger (see Fig. 3, dashed line in the bottom panel). The second effect is the radial 
variation of the efficiency of the detection of Subs galaxies. Very simple modeling shows 
that the efficiency of the selection method decreases towards the cluster center. As a 
result, both the number density profile and the vdp of the Subs class are affected. A 
simple first-order correction of these profiles leads to a /^'-profile which no longer implies 
tangential orbits (see Fig. 3, dash-dotted line in the bottom panel). A better under- 
standing of the orbital characteristics of subclusters thus requires either an improved, 
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radially-unbiased, selection method of galaxies in subclusters, or at least an improved 
modeling of the biases. 

4. Summary 

The Early-class galaxies (ellipticals and SO galaxies) in clusters have a distribution 
of line-of-sight velocities that is consistent with an isotropic orbit distribution. These 
galaxies were therefore used to derive the mass profile of the ENACS ensemble cluster. 
Using this mass profile, we constrained the orbits of cluster galaxies of other classes. From 
the distribution and velocity dispersion profile of the S-early-class we conclude that these 
may also have isotropic orbits, although they show an apparent radial anisotropy around 
w 0.45 r2oo- The latter may be a result of the process by which the majority of the S-early- 
galaxies are thought to have transformed into SO galaxies through impulsive encounters. 
The fact that the local density around S-early-galaxies is smaller than that around SO 
galaxies provides further support for this interpretation. The S-late-class has mildly radial 
orbits outside ~ 0.7r 2 oo, and are probably still falling in. Their absence from the central 
region most likely indicates that they are tidally disrupted on their crossing of the cluster 
core. Finally, the galaxies in substructures apparently are on tangential orbits, although 
isotropic (or mildly radial) orbits cannot be excluded at present. 
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